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Cn , The recent developments on the computation of non-resonant corrections to the e + e — > 

i—i ■ W + W~bb cross section in the top-antitop resonance region are reviewed in this talk. 

Non-resonant production of the final state W + W~bb starts to contribute at NLO in 
the nonrelativistic power-counting v ~ a s ~ y/awN- The corrections induced by non- 
resonant effects reduce the cross section in the top-antitop resonance peak region at a 
level which is comparable to the expected experimental precision at the future linear 
collider, and are thus relevant for a high-precision top mass and width determination. 
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^' 1 Introduction 

P-i- 

D ■ A next generation e + e linear collider such as the International Linear Collider (ILC) oper- 

tG [ ating at energies around the top-antitop threshold will allow us to measure the top-quark's 

properties with unprecedented precision. The top-quark mass is currently known from direct 
production at the Fermilab Tevatron (and soon at the Large Hadron Collider) with a preci- 
>> i sion > 1 GeV. From a threshold scan of the e + e~ — > tt cross section at the ILC, however, an 

order of magnitude improvement in the precision can be achieved experimentally [Tj . Aside 
from determining a fundamental parameter of the Standard Model, its precise knowledge is 

f->) | important for precision tests of the Standard Model and its extensions. Other characteris- 

tics of the top quark such as its width and Yukawa coupling provide information about its 
coupling to other particles and the mechanism of electrowcak symmetry breaking. For these 

r»o . reasons top-quark pair production near threshold in e + e~ annihilation has been thoroughly 

investigated following the nonrelativistic QCD (NRQCD) approach, which treats the lead- 
ing colour-Coulomb force between the nearly on-shell top and antitop exactly to all orders 
in perturbation theory. In this framework, where the strong coupling a s is of the same 
order as v, the small relative velocity of the top and antitop, QCD corrections to the total 
cross section are known at the next-to-next-to- next-to-leading order (NNNLO) [2HI]j and 
higher-order logarithms (ev s logv) fe have been resumed to next-to-next-to-leading logarith- 
mic order (NNLL) [5|. The NNLL prediction, which had a normalization uncertainty in the 
total cross section > 6%, has been recently updated jB], showing a substantial reduction of 
the renormalization scale dependence. An analysis collecting all the pieces that contribute 
to the full NNNLO total cross section is still pending (see [3] for the latest summary), which 
would provide information about the convergence of the non-logarithmic piece of the QCD 
corrections. 

Here we focus on sublcading corrections of electroweak origin. The top quark is unstable 
with a significant width Ft of about 1.5 GeV due to the electroweak interaction. Once 
the top width is included, due to top decay, the physical final state is W + W~bb - at 
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least if we neglect the decay of top into strange and down quarks, as justified by Vtb ~ 1, 
and consider W bosons as stable. The W + W~bb final state can also be produced non- 
resonantly, i.e. through processes which do not involve a nearly on-shell tt pair. The 
latter effects are not included in the standard nonrclativistic treatment used to compute the 
dominant QCD corrections, which can only describe the resonant process with a nearly on- 
shell tt intermediate state, but can be accommodated using the effective field theory (EFT) 
formalism to describe pair production of unstable particles at threshold [7HS]. Adopting a 
counting scheme where oew ~ o%, the leading non-resonant and off-shell effects are NLO 
for the total cross section, since there is an additional power of qew but no phase space 
suppression, hence the relative correction is cvew/v ~ cn 8 . At NLO they can be classified as 
part of the electroweak corrections to the e + e _ — > W + W~bb cross section, and are dominant 
compared to the purely resonant electroweak effects, which first contribute at NNLO [5]. The 
identification of non-resonant production of the final state W + W~bb as a pure electroweak 
effect does not hold anymore at NNLO, where gluon corrections to the non- resonant part 
have to be considered. 

In this talk I discuss the status of the calculation of the non-resonant contributions to 
the e + e~ — > W + W~bb process in the tt resonance region. The results at NLO are fully 
known for the total cross section as well as including invariant-mass cuts on the Wb pairs. 
At NNLO, interesting conceptual issues regarding the interplay between resonant and non- 
resonant parts of the calculation arise. These are discussed in Sec. [3] The computation of 
the full set of NNLO non-resonant corrections represents a much more difficult task that has 
not yet been attempted. An alternative approach, named "NRQCD phase space matching" , 
to compute the non-resonant effects beyond NLO entirely through calculations in NRQCD 
has been proposed |10H11) . which works if moderate invariant-mass cuts on the Wb pairs 
are applied. The concepts of the phase space matching and the results from this method 
are reviewed in the last section. 

2 Non-resonant NLO electroweak contributions 

The cross section for the e + e~ — > W + W~bb process is obtained from the W + bW~b cuts of 
the e + e~ forward-scattering amplitude. In the energy region yj~s ps 2mt the amplitude is 
dominated by the production of resonant top quarks with small virtuality. A separation of 
resonant and non-resonant effects can be consistently achieved within the unstable-particle 
effective field theory [7]. This allows us to integrate out hard modes (~ mt) and represent 
the forward-scattering amplitude as the sum of two terms, 

iA = Y. C p k)C p } f d i x(e-e + \T[iO^(0)iO^(x)]\e~e + ) 

t ) (e-e + \ i o£ ) (0)\e-e+). (I) 
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The matrix elements in (Q]) are evaluated in the "low-energy" effective theory, which 
includes elements of soft-collinear and nonrelativistic effective theory. The first term on the 
right-hand side of ([!]) describes the production of a resonant tt pair in terms of production 
(decay) operators O p (x) (O p (x)) with short-distance coefficients C p . The second term 
accounts for the remaining non-resonant contributions, which in the effective theory are 
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Figure 1: Two-loop forward-scattering amplitude diagrams with tbW + cuts. tbW cuts and 
symmetric diagrams are not shown. 

described by four-electron production-decay operators Ae ' . The coefficients C±J originate 
from the hard contributions of the e + e~ forward-scattering amplitude. The hard momentum 
region expansion dictates that the top-quark self-energy insertions are treated pcrturbatively, 
since the top lines are formally far off-shell, p 2 — m 2 ~ 0(m 2 ) 3> S(pf). Accordingly 

(k) 

the calculation of the coefficients C\ e is performed in fixed-order perturbation theory in 
the full electroweak theory with no resummation of self-energy insertions in the top-quark 
propagator [5] and supplemented with an expansion of the amplitudes near threshold (in 
8 = s/(4m 2 ) — 1 ~ v 2 ). The leading imaginary parts of C ie arise from the cut two-loop 
diagrams of order a^ w shown in Fig. [T] The corresponding contribution to the cross section 
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Technically, this simply amounts to the calculation of the spin-averaged tree-level processes 
e + e~ — > tW~b and e + e~ — > tW + b with no width supplied to the intermediate top-quark 
propagators. Instead, the divergence from the top-quark propagators going on-shell is regu- 
larized dimensionally. Details on the computation and integral representations of the result 
for ([2]) can be found in [12]. The results of [12] for the NLO non- resonant contribution 
to the total cross section have been confirmed recently by an independent calculation |13j . 
The latter is based on an expansion in the parameter p — \ — M\y / 'mt rj 0.53. which for 
individual diagrams requires to consider several orders in p or to use Pade approximants to 
reach a precise numerical agreement with the integral representation of |12j . For the sum of 
all diagrams, however, the leading order in p gives an approximation which differs from the 
exact result by less than 5%. 

Through the computation of the four-electron matching coefficients loose cuts (~ m t ) 
on the bW + and bW~ invariant masses can be incorporated easily, as it has been discussed 
in the context of VF-pair production near threshold Q3] . The result obtained in fT2] covers 
the case of symmetric cuts on the invariant mass of the bW subsystems {p 2 w ) of the form 



m t - AM t < ^Jp^w < m * + AM t , 
for AMt ^> Tt, up to the total cross section (AM( iinax = m t — Mw)- 



(3) 
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Figure 2: Relative sizes of the QED and non-resonant corrections with respect the tt LO 
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cross section in percent: ctqed/ ** (upper solid blue line), a 



(0) 



(i) 



noil- res 



/a (0) 
/ u tt 



for the total cross 

section (lower solid red line) and AM t = 15 GeV (lower dashed red line). The relative size 
of the sum of the QED and non-resonant corrections is represented by the middle lines, for 
AM t>max (solid) and AM t = 15 GeV (dashed). We have chosen a s (30GcV) = 0.142 and 
m t = 172 GeV. 



2.1 Results 

The plot in Fig. [2] displays the relative sizes of the NLO electroweak corrections with respect 
to the LO result for the e + e~ — » W + W~bb cross section, which includes the summation 
of Coulomb corrections. The QED contribution represents a correction of about 2% above 
threshold and rises to a maximum of 7% just below the peak, while the non-resonant con- 
tributions give a constant negative shift of about 3% above threshold. Below threshold the 
relative size of the non-resonant corrections is very large, since the LO result rapidly van- 
ishes, reaching up to 19%. Hence below threshold they represent the leading electroweak 
correction to the total tt cross section. We observe a partial cancellation of the QED and 
non-resonant corrections in the peak region and at energies above. A sensitivity to the 
invariant-mass cut AM t in the bW + and bW~ subsystem enters first at NLO through the 
non-resonant contributions. Restricting the available phase space for the final-state particles 
by tightening the invariant-mass cuts AM t makes the non-resonant contributions even more 
important. This is shown by the dashed lines in Fig. O corresponding to AM t = 15 GeV. 
The non-resonant correction amounts to a negative shift of 27-35 fb for yG in the interval 
(338,350) GeV. 



3 Non-resonant effects beyond NLO 

3.1 Finite-width divergences 

There is an interesting conceptual issue concerning the resonant part of the QCD calculation 
of the tt cross section. At NNLO it exhibits an uncanceled ultraviolet divergence (here 
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Figure 3: Tree-level double-resonant diagrams in the full and effective theories. Double lines 
are used for representing top quarks. The (red) dashed line denotes that we extract the 
imaginary part of the forward scattering amplitude or, equivalently, that we perform the 
phase space integration over the particles in the cut. 



regulated dimcnsionally) 



a,T 



u t i ex oc , (4) 

which arises from the logarithmic overall divergence in the two-loop nonrelativistic corre- 
lation function, whose imaginary part gives the cross section. The overall divergence is 
polynomial in the nonrelativistic energy E of the top quarks, but contributes to the cross 
section if T t ^ 0, since the correlation function is evaluated at complex values E — >• E + iT t . 
It can be shown that these UV-divergences originate because the unstable particle prop- 
agators describing the top and antitop in the EFT allow for contributions to the forward 
scattering amplitude from intermediate states which have arbitrarily large invariant masses. 
The inclusion of a finite width in the top propagator changes the high-energy behaviour 
of its imaginary part and makes the phase space integration extend to infinity. To illustrate 
the latter consider the Born cross section diagram in the EFT, Fig. [3l The cutting rules for 
the top lines imply extracting the real part of the NRQCD propagator: in the stable case 
(iT t /2 — > ie) we have 
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which imposes a fixed dispersion relation p° = p 2 /2mt, while the real part of the unstable 
propagator, 

r t /2 
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yields a Brcit-Wigner distribution with support in the entire (p°, |p|) plane. 

From the point of view of the full theory, already taking the nonrelativistic limit makes 
the phase space integration extend to infinity. Consider the full-theory double-resonant 
diagram of Fig. [3] with center-of-mass momentum q = (2m t + E, 0). When the quark lines 
are close to their mass shell the 4-particle phase space integration reduces to 



&tt 
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\m t E + 2m t po p 2 + ie\ \m t E - 2m t po - p 2 



(7) 



where in the denominator we have retained the leading order term in the nonrelativistic 
expansion of the top and antitop off-shellness, (p 2 f — m 2 ), with p 2 tl = (q/2 ±p) 2 . The 
r 2 factor in the numerator arises from the phase space integration of the bW + and bW~ 
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subsystems because in the nonrelativistic limit the top and antitop quark are effectively on- 
shell. We notice that the integration limits in Eq. ([7}, which in the full theory computation 
are cut off by the top mass, become infinite in the nonrelativistic limit. The boundaries of the 
phase space are determined by step functions in the phase space measure, and the arguments 
of these functions also have to be expanded according to the nonrelativistic power-counting. 
In the limit m t ^> (p°, |p|) the step functions do not depend on (p° , |p|) any longer, and 
are satisfied trivially, thus allowing for an infinite integration region in these variables. The 
corresponding NRQCD amplitude (diagram on the right in Fig. [3]) reproduces Eq. (JT)) when 
r t is treated as an insertion. The NRQCD power-counting, however, tells that the r t ~ E 
term needs to be rcsummed as part of the top propagator, thus effectively replacing ie by a 
term proportional to iT t as in Eq. ((6]). 

Despite the integration limits, the tree-level integration in Eq. (JT]) is finite. However the 
integrand becomes more sensitive to large momentum regions once we include relativistic 
corrections ~ p 2 /m 2 . Using dimensional rcgularization these sublcading contributions can 
lead to 1/e singularities if the high energy behaviour of the EFT phase space integration is 
logarithmic. The lesson from this is that the pure resonant result alone that is usually shown 
in the literature is inconsistent theoretically and must be supplemented with additional 
short-distance information from the systematic calculation of the e + e~ — > W + W~bb process. 

3.2 Towards an evaluation of the NNLO non-resonant contributions 

Since the full-theory calculation is finite, the UV divergence in the NNLO resonant part must 
cancel with an infrared divergence that appears in the non-resonant term in unstable-particle 
effective theory from diagrams corresponding to off-shell top-quark decay, as discussed in [4] . 
The NNLO contributions to er non _ rcs are given by the 0(a s ) corrections to the diagrams 
/ii-io in Fig. Q] (an example is provided by Fig. [4j. Parts of the latter consist of gluon radi- 
ation diagrams, which thus contribute first at NNLO in the non-resonant part, in contrast 
to what happens in the resonant part, where ultrasoft gluon radiation is a N 3 LO effect. 

Divergences in the non-resonant part arise when 
the top (or antitop) propagators go on-shell [12] . 
This is a consequence of the hard momentum re- 
gion expansion, which sets T t = in the top-antitop 
propagators and forces to retain the leading order 
term of the upper kinematic limit in p 2 = p\ w '- 

Pt.max = m \ +0(S). Integrating over all other kine- Figure 4: NNLO contribution to 
matic variables but p 2 , the non-resonant contribu- 
tions involve integrals of the form 
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where the endpoint singularity at p 2 = m 2 is regularized in d = 4 — 2e dimensions. At NLO, 
only the diagram h\ has an endpoint divergence, with n = 3/2, and the result is therefore 
finite in the limit e — > 0. Gluon corrections, on the other hand, yield additional half-integer 
powers of (m 2 — p 2 ) and, in particular, n = 1 contributions which will generate 1/e terms 
of the form ((4|). For example, the gluon vertex correction to diagram hi shown in Fig. |4] 
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yields non-zero contributions with n = 2,|,l,i,..., and one can check that the resulting 
divergences cancel those in the resonant part of the full-theory diagram. Clearly, the com- 
plexity associated with the evaluation of the full NNLO set of non-resonant corrections calls 
for a numerical implementation using Monte Carlo methods for the phase space integration. 
Such an implementation requires that the end-point divergences occurring in the diagrams 
are identified and subtracted from the amplitude, together with the soft-collinear divergences 
that show up at NNLO due to gluon radiation. The complete set of NNLO non-resonant 
singularities shall be available soon (15] . 

It should be mentioned that the leading term of the expansion in p = 1 — Myy /m t of the 
NNLO non-resonant contribution to the total cross section has been given in Ref . [13] . The 
accuracy of such an approximation, though, cannot be estimated reliably without further 
knowledge of the size of subleading terms in the p-expansion. 



4 Phase space matching 

An alternative approach to account for non-resonant effects has been developed in paral- 
lel [IQlIll] that includes the effects of invariant-mass cuts on the Wb pairs entirely through 
calculations in an extended version of NRQCD which accounts for unstable particle effects 
systematically [H] . In this approach, the UV-divergences from the resonant matrix elements 
are compensated by imaginary counter-terms associated with the (e + e~)(e + e~) forward- 
scattering operators, which thus acquire an imaginary anomalous dimension and sum large 
logarithms of the top velocity In Ref. [IT] we demonstrate that for moderate top invari- 
ant mass cuts of the form ([3j with AM t ~ 15 — 35 GeV, the matching conditions of the 
(e + e~)(e + e~) forward scattering operators are dominated by the NRQCD phase space con- 
tributions, i.e. they can be computed from the difference between the (potentially) divergent 
NRQCD phase space integrations without any cuts and the ones with the cuts in Eq. ([3]) 
being imposed. This is because using the MS scheme in NRQCD diagrams involving the 
unstable top propagator of Eq. ([6]) largely overestimates the contributions from unphysi- 
cal phase space regions that are paramctrically away from the potential, soft and ultrasoft 
regions that can be described by NRQCD. Thus the main numerical effect of the "phase 
space matching" procedure is obtained by removing these unphysical contributions and can 
be carried out within NRQCD itself. 

An important conceptual aspect of the invariant mass cuts defined in Eq. ([3]) is that 
already for moderate cuts AM t ~ 15 — 35 GeV the cut Ap on the nonrelativistic (anti)top 
three-momentum p is Ap ~ \/2m t AM t ~ 100 GeV, and thus represents a hard scale of 
the order m t . This justifies the implementation of the phase space effects into the matching 
conditions of the Wilson coefficients. Since A = ^/2m t AMt is paramctrically of order vtit 
we use for our bookkeeping the counting A ~ mt. In this counting scheme the phase space 
constraints are incorporated through the NRQCD Wilson coefficients. On the other hand, 
numerically the scales AMt and A are sufficiently below the top mass scale such that all tt 
phase space configurations that pass the invariant mass constraint can still be adequately 
described by NRQCD. This fact is crucial for the phase space matching method briefly 
outlined in the following. 

Consider first the case without QCD effects (a s = 0) . For the determination of the Wilson 
coefficients CyK we need to know the result for the inclusive cross section Cj"*^ (A) with 
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invariant mass constraints q In the common approach to matching computations, (Jj"*J" (A) 
is computed in the full relativistic theory. After the result is expanded nonrelativistically 
using the counting v ~ y'ciiw, one can identify the pieces belonging to the Wilson coefhcients 

Cy/ A - On the other hand, as mentioned above, the tt phase space regions passing the 
invariant mass cuts can be determined within the nonrelativistic expansion. We therefore 
write the expression for the inclusive cross section as a sum of two terms, 

C c r°(A) = a-=° CD (A) + O m =°(A) • (9) 

Here, ct^j^q ) cd is the cross section computed from NRQCD Feynman rules with the (anti)top 
invariant mass constraints being applied for the phase space integration. The parameter A 
is related to the invariant mass cut AMt and we use the formal counting A ~ mt according 
to the discussion above. In this computation the (e + e~)(e + e~) forward scattering operators 
do not contribute, and the resulting expressions are just the nonrelativistic expansions of 
full theory squared matrix elements containing the square of the double resonant diagram 
e + e~ — >• tt — >• W + W~bb (see Fig. EK) an d the interference of the double resonant diagram 
with the diagrams for e + e~ — > W + W~bb having only either the top or the antitop in 
intermediate stages (see Fig. [5}d and c for typical diagrams). The contributions to the 
Wilson coefhcients Cy, A that result from cr^|^ CD (A) are local (i.e. energy-independent) 
and have the form 



- y 

A Z-> 

n,k=0 



m t T t \ " ( A 



(10) 



where the first term in the expansion, proportional to r t /A, arises from the leading order 
diagram for the phase space matching computation (right diagram in Fig. [3J and gives 
the dominant phase space correction to the tt cross section. In the counting A ~ m tl it 
constitutes an 0(v 2 ) correction (NLL0), since T t /m t ~ v 2 . The A 2 /m 2 terms arise from 
insertions of operators that are higher order in the nonrelativistic expansion, and since they 
are formally of order unity, can lead to power-counting breaking contributions. However, 
we find that the numerical effects of the power-counting breaking contributions are very 
small and do not spoil the nonrelativistic expansion. This is partly due to the fact that the 
phase space cutoff A is sufficiently smaller than the convergence radius of the nonrelativistic 
expansion. The phase space matching procedure can thus be implemented for values of the 
invariant mass cut in the range m t r t C A 2 < to 2 , but not for the total cross section where 
A is numerically of order m* . 

The remainder contribution of the inclusive cross section, v"^ (A) accounts for all other 
contributions to the full theory matrix element. This includes for example pure background 
e + e~ — > W + W~bb diagrams, see Fig. [5jl for a typical diagram, and also the square of the 
single-top diagrams in Figs. [5}d and c. One can check that the first two terms in the phase 
space matching series (fTU)) agree with the two first terms in the expansion in A 2 /m 2 of the 
full NLO non-resonant result. This means that the remainder contributions, not calculable 



a The Wilson coefficients Cy,. correspond to the C\J coefficients of Eq. {U, following the notation of |llj . 
which is explained in their Sec. II. 

b In the NRQCD approach were logarithms are summed systematically the counting LO, NLO. . . is 
replaced by leading logarithmic (LL) order, next-to-leading logarithmic (NNL) order and so on. 
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-> tt — > W + W~bb. (b,c) Typical single-resonant 
full theory diagram for e + e~ —} W~tb — > W + W~bb. (d) Typical full theory diagram for 
W + W~bb without top or antitop quarks as intermediate states. 



Figure 5: (a) Full theory diagram for e + e 

o + 



within the phase space matching approach, first contribute with r t A 3 /ra^ terms, and should 
therefore be small for the range of A used in |llj . Indeed, the remainder contribution was 
determined in [TT] by computing the full e + er — »• W + W~bb cross section at tree-level with 
MadGraph, and shown to be smaller than 5 fb; it can thus be neglected in view of the 
expected experimental precision expected at a future linear collider. Since no kinematic or 
dynamical enhancement is expected for the QCD corrections to the remainder contributions, 
we assume these are also negligible when a s ^ 0. This simplifies the computations substan- 
tially and makes the determination of higher order QCD corrections feasible in the phase 
space matching approach. A calculation of the full NNLO non-resonant contributions along 
the lines of Sec. 13.21 shall confirm this assumption. 

The QCD corrections introduce powers of (a s m t /A) in the phase space matching series 
([TO]) . As far as the a s -cxpansion is concerned, we have found that the N 3 LL (0(a 2 s )) 
corrections to the phase space matching contributions need to be determined to meet the 
experimental precision expected at a future linear collider. An important part of these N 3 LL 
corrections have been computed analytically in [TT] , together with the complete set of NLL 
and NNLL ones. It is interesting to note that at N 3 LL order one has to include also the phase 
space matching for (e + e~)(tt) top production operators. This is because at higher orders 
in the loop expansion, one has to account for the phase space matching contributions of 
subdiagrams to remove non-analytic matrix element terms from the matching equations and 
to achieve that the matching coefficients are analytic in the external energy. The procedure 
of carrying out the phase space matching is thus analogous to the common matching and 
rcnormalization methods for stable particle theories. 



4.1 Results 

In Fig. |5] we show the phase space matching corrections to the inclusive tt threshold cross 
section up to N 3 LL order, for invariant mass cuts of AM t = 15 GeV and AM t = 35 GeV, 
and compare them with the rest of electroweak corrections. The QED effects arise from the 
electromagnetic correction to the QCD Coulomb potential and the one-loop QED match- 
ing correction to the Wilson coefficient of the tt current. The hard one-loop electroweak 
corrections have been obtained in Ref. |16j . The type-1 finite lifetime corrections represent 
all finite lifetime corrections which are not related to phase space constraints. They consist 
of the corrections generated by the imaginary interference matching coefficient, the time 
dilation corrections to the Green function, both known at NNLL order, and the contribu- 
tions from the renormalization group summation of the NLL phase space logarithms [5]. 
The QED, hard electroweak and type-1 finite lifetime corrections do not depend on phase 



LCWS11 



b 




342 344 346 

V7 (GeV) 



342 344 346 

V7 (GeV) 



Figure 6: Sizes of the different contributions to the inclusive cross section arising from 
electroweak interactions as a function of the total c.o.m energy for AM t = 15 GeV (left) and 
AM t = 35 GeV (right): (green line) NNLL hard one- loop electroweak effects from Ref. [T6] . 
(red line) NNLL finite lifetime corrections from Ref. |5], (blue line) NNLL QED effects, 
and phase space matching corrections at NLL, NNLL and N 3 LL (dotted, dashed and solid 
black lines, respectively). The blue dotted lines correspond to the expected experimental 
uncertainties at the LC. The values used for input parameters are those of Ref. [llj . 



space restrictions and are therefore identical in both panels. In Fig. |6]the blue dotted lines 
represent a rough estimation of the expected experimental uncertainties at a future linear 
collider consisting of an energy-independent error of 5 fb and a 2% relative uncertainty with 
respect to the full prediction, both being added quadratically. 

We see that the QED (blue lines) and the type-1 finite lifetime corrections (red lines) are 
sizeable (at the level of 40 fb) only in the peak region just below y/s = 2m t . Due to their 
different signs the QED corrections and the type-1 finite lifetime corrections cancel each 
other to a large extent in the peak region. The hard electroweak corrections (green lines) 
represent a multiplicative factor of -1.2% to the total cross section and are therefore very 
small below the peak and at the level of 12-13 fb above the peak region. We see that the 
phase space matching contributions represent the largest of the four classes of electroweak 
effects. In contrast to the other classes of electroweak effects they do not decrease strongly 
for energies below the peak region. For AMt = 15 GeV the N 3 LL phase space matching 
contributions amount between —85 and —65 fb and for AM t = 35 GeV they are between 
—45 and —35 fb. The overall size of the phase space matching corrections decreases for larger 
values of the top invariant mass cut AM t . We emphasize, however, that the results obtained 
within the phase space matching approach are valid only for moderate values of AM t in the 
region between 15 and 35 GeV. For invariant mass cuts below 15 GeV the phase space 
constraints arc not related anymore to hard effects and for invariant mass cuts substantially 
above 35 GeV matching contributions that need to be computed from full theory diagrams 
have to be included. The relatively flat behavior of the phase space matching contributions 
is related to the fact that the dominant phase space matching contributions are energy- 
independent. The small linear dependence on y/s is related to the y/s dependence of the 
virtual 7 and Z propagators of the basic e + e~ — > tt process and the peak-like structure comes 
from an imaginary phase space matching contribution to the (e + e~)(tt) top pair production 
operator which enters the N 3 LL inclusive cross section in terms of a time-ordered product. 
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Figure 7: Total inclusive top pair production cross section from NRQCD: starting from the 
pure QCD NNLL prediction (black dashed line), we add step- by-step the QED corrections 
(blue line), the hard electroweak corrections (green line), the type-1 finite lifetime corrections 
(red line) and the N 3 LL phase space corrections (black solid line) for AM t = 35 GeV. 



The results for the NLL, NNLL and N 3 LL phase space marching contributions displayed in 
Fig. [6] show that the expansion related to the phase space matching procedure is particularly 
good for larger values of AM t and still well under control for AAI t = 15 GeV. We note that 
the rather small size of the NNLL corrections (difference of black dotted and dashed lines) for 
AM t = 15 GeV arises from a cancellation between different independent NNLL corrections. 
In Fig. [7] the size of the four different types of electroweak corrections is shown for 
predictions of the total inclusive cross section for AM t = 35 GeV. We again see that the 
phase space matching contributions exceed by far the other electroweak corrections. The 
phase space matching contributions are between —85 and —35 fb for invariant mass cuts 
AM t between 15 and 35 GeV and are essential for realistic theoretical predictions. In the 
peak and the continuum region (t/s > 2mt) they amount to 6 to 10%. They are particularly 
important in the region below the peak (t/s < 2mt) where the cross section decreases and 
the unphysical off-shell contributions of the NRQCD tt phase space become dominant. Here 
the phase space matching contributions can amount to more than 50%, and they ensure 
that the cross section has the correct physical behavior. 
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